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Abstract 

Phenomenological orbital-polarization (OP) terms have been repeat- 
edly introduced in the single-particle equations of spin-density-functional 
theory, in order to improve the description of orbital magnetic moments 
in systems containing transition metal ions. Here we show that these 
ad hoc corrections can be interpreted as approximations to the exchange- 
correlation vector potential A^c of current-density- functional theory (CDFT). 
This connection provides additional information on both approaches: Phe- 
nomenological OP terms are connected to first- principles theory, leading 
to a rationale for their empirical success and a reassessment of their limi- 
tations and the approximations made in their derivation. Conversely, the 
connection of OP terms with CDFT leads to a set of simple approxima- 
tions to the CDFT potential A.xc, with a number of desirable features 
that are absent from electron-gas-based functional. 
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1 Introduction 



Magnetic phenomena arising as a consequence of the spin degrees of freedom and 
the antisymmetrization of the wave function are ubiquitous in quantum chem- 
istry and physics, and are routinely treated within both density-based and wave- 
function-based approaches. Magnetic phenomena arising from orbital currents, 
on the other hand, are not automatically accounted for by standard methods 
for electronic-structure calculations. The present paper is dedicated to explor- 
ing links between a phenomenological and a first-principles approach to orbital 
magnetism in density- functional theory (DFT). These links shed new light on 
both approaches, and pave the way for new applications of each. 

In Sec. [21 of the present paper, we briefly recall a number of distinct ways 
in which orbital magnetism can occur in a many-electron system. In Sec. |31 
we describe two computational approaches to orbital magnetism. Section 13.11 
is devoted to current-density functional theory (CDFT), as a first-principles 
method for the description of orbital magnetism in DFT [1,2]. Independently 
of CDFT, a variety of phenomenological orbital polarization (OP) terms were 
proposed as add-ons to the single-particle equations of spin-density- functional 
theory (SDFT) [3-7]. These terms are quite commonly used in band-structure 
calculations of magnetic solids, but are much less known in quantum chemistry. 
We devote Sec. l3.2l to a quick description of the motivation and form of several 
such add-on terms. 

In Sec. 21 we show that these phenomenological terms can be related to the 
first-principles approach, both for single-particle energies (Sec. 14. ip and total 
energies fSec. 14. 2j ). This connection benefits both sides: Phenomenological OP 
terms are connected to first-principles theory, leading to a rationale for their 
empirical success and a reassessment of their limitations and the approximations 
made in their derivation. Conversely, the connection of orbital-polarization 
terms with CDFT leads to a set of simple approximations to the CDFT vector 
potential A^c of systems containing open-shell atoms, with a number of desirable 
features that are absent from electron-gas-based functionals. 

Section [5] contains our conclusions. 

2 Orbital magnetism in many-electron systems 

There are at least four conceptually distinct ways in which orbital magnetism 
can appear in a physical system. One is the presence of external magnetic fields 
B(r), whose vector potential A(r) enters the Hamiltonian via the usual minimal 
substitution in the kinetic energy 



A second way in which orbital magnetism can appear is due to current- 
current interactions, which are part of the Breit interaction [8,9] and therefore 




(1) 
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a relativistic effect. The nonretarded part of this interaction is 

-ijo''-l ^r' '■'f^ = -f / "■'"j-.i') ^ A„ W. (2) 

which describes the Hartree-hke couphng of currents to the self-induced vector 
potential, corresponding to the Amperian currents of classical electrodynamics. 

Third, in the presence of spin-orbit coupling a nonzero spin magnetiza- 
tion can induce an orbital magnetization. Orbital magnetic moments induced 
by spin-orbit coupling can be treated within relativistic DFT, and are often 
described within SDFT by adding a spin-orbit coupling term to the Hamil- 
tonian. Such magnetic moments become important, e.g., in magnetic solids, 
where spin-orbit coupling produces phenomena such as magneto-crystalline an- 
isotropy [10, 11] or magneto-optical effects such as dichroism [12, 13]. 

Finally, and most intriguingly. orbital magnetism can also occur sponta- 
neously in a system with pure Coulomb interactions, if the system minimizes 
its energy in a current- carrying state [2, 14]. The resulting currents are in prin- 
ciple functionals of the charge density, because the original formulation of the 
Hohenberg-Kohn theorem applies in the absence of magnetic fields, but (S)DFT 
provides no explicit prescription how to calculate the spontaneous orbital cur- 
rents and their effect on observables. This situation has changed with the advent 
of nonrelativistic current-density-functional theory, developed by Vignale and 
Rasolt [1,2], which describes spontaneous currents by introducing in the Kohn- 
Sham equations a self consistent exchange-correlation (xc) vector potential Axc^ 
which can be nonzero also in the absence of external magnetic fields and of rel- 
ativistic effects. 

In addition to possibly appearing spontaneously, A-xc also becomes nonzero 
as soon as currents are induced by one of the other three mechanisms. In this 
case, it constitutes a correction to the external or internal vector potentials, or 
spin-orbit terms, already present in the Hamiltonian. 

CDFT has been applied to the calculation of the effects of orbital magnetism 
in atoms [15-17], quantum dots [18-20], molecules [21,22] and solids [14,23-25]. 
In these applications it has become clear that the main challenge of CDFT at 
the present stage of its development is the construction of reliable and computa- 
tionally viable approximations to the exchange-correlation vector potential A.xc- 
The present paper explores the possibility to develop such approximations by 
drawing on analogies to independently developed phenomenological approaches 
to orbital magnetism. 

3 Phenomenological and first-principles treat- 
ment of orbital magnetism in density-functional 
theory 

This section collects, without proofs (which can be found in the original lit- 
erature), the key equations of both the first-principles (CDFT) and the phe- 
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nomenological (OP) approach to orbital magnetism in DFT. As our aim is to 
unravel a connection between the two, we limit ourselves to providing a brief 
overview, highlighting the aspects that will allow us to identify the phenomeno- 
logical approach as a well-defined approximation to the first-principles one. 



3.1 Current-density-functional theory 

To briefly describe the formalism of nonrelativistic CDFT, we first recall the 
form of the traditional Kohn-Sham (KS) equation of DFT, 



Here the effective single-particle potential, Vs{r), is defined as 



Vs{r) = u(r) + vnir) + Vxdr), 



(3) 



(4) 



where v{r) is the external potential, t'H(r) the Hartree electrostatic potential, 
and Vxci^) the exchange-correlation potential, in which the entire complexity of 
the many-body problem is hidden [26,27]. 

By construction, the single-particle orbitals 4'k solving the eigenvalue prob- 
lem ^ reproduce the density of the interacting system via 



(5) 



On the other hand, the paramagnetic current density, 



2rm ^ — ' 



(6) 



and the orbital magnetic moment 



KS 



^ I d^r(zxr).jf«(r) 



(7) 



following from these orbitals, are, a priori, not guaranteed to have any relation 
with the true current density and magnetic moment of the interacting system. 

The corresponding equations of CDFT have a slightly more complicated 
form: 



2m ' " " ' 



where an upper index 'c' denotes CDFT, 



(8) 



(A(r)2-A,(r)2 



2mc^ 

D^(r) v{r) + VH{r)+v''^dr), 



(9) 
(10) 



5 



and 



A,(r) = A(r) + A^e(r). 



(11) 



Here v^^ and A^c are the exchange-correlation scalar and vector potentials of 
CDFT, respectively [1,2]. A.xc, in particular, is a gauge invariant functional 
of the densities n(r) and jp(r), written Axc[n,ip]{r). By setting A^c = 0, one 
recovers from CDFT the equations of (S)DFT in an external vector potential. 
The novel feature of CDFT is that A^c accounts for the orbital degrees-of- 
freedom in the single-particle equations, even in the absence of external fields, 
and allows calculation of orbital currents directly from a set of Kohn-Sham-type 
equations. 

The single-particle orbitals tpk solving the more complicated eigenvalue prob- 
lem ^ reproduce by construction both the density 

n{r)^Y.^l{r)Mr), (12) 

k 

and the paramagnetic current density 

= [^^(^)^^'»(^) " (V^:(r))^,(r)] (13) 

k 

of the interacting many-body system. The correct orbital magnetic moment can 
then be obtained from 

(4)-^ I dV(zxr).j,(r). (14) 

The paramagnetic current alone is not gauge invariant, but the gauge invari- 
ant orbital current, jorb(r), is simply obtained from 

jorb(r)=jp(r)-^n(r)A(r). (15) 

mc 

An important property of the effective potential of DFT, Vs(j), which is not 
shared by the CDFT potentials V/(r) and As{r) [28], is its uniqueness: for any 
given ground-state density n(r) there is, up to an irrelevant additive constant, 
at most one such local multiplicative potential [26] . 

Any CDFT calculation requires an approximation for the current dependence 
of the Exc functional. For the homogeneous three-dimensional electron liquid in 
strong uniform magnetic fields, the exchange energy is known exactly [29], and 
the correlation energy has been calculated within the random-phase approxima- 
tion [30] and the self-consistent local-field corrected Singwi-Tosi-Land-Sjolander 
scheme [31]. In weak fields, where linear-response theory applies, the exchange- 
correlation energy can be expressed in terms of the magnetic susceptibility, for 
which many-body calculations are available from [32]. The exchange-correlation 
energy of two-dimensional electron liquids in uniform magnetic fields has been 
much studied in the context of the fractional Quantum Hall effect in quasi-two- 
dimensional semiconductor heterostructures [33] . 
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3.2 Orbital-polarization terms 



In an independent take on the problem of orbital magnetism, originating with 
Brooks and collaborators [3,4], phenomenological orbital-polarization terms are 
introduced in the formal framework of SDFT. This effort is motivated by ob- 
serving that SDFT does account for Hund's first rule (dealing with spin angular 
momentum) and can be extended to account for Hund's third rule (due to spin- 
orbit coupling), but does not contain an obvious ingredient corresponding to 
Hund's second rule (dealing with orbital angular momentum). 

In the absence of external magnetic fields, the exact exchange-correlation 
functional of SDFT doubtlessly would predict the correct orbital magnetic mo- 
ments and orbital currents, as functionals of the spin and charge densities, but 
these functionals are not known. In practice, such quantities are therefore of- 
ten calculated directly from the orbitals of SDFT. We note that there are two 
separate issues here: One is the use of SDFT orbitals to calculate orbital mag- 
netic moments, although these orbitals are constructed to reproduce only the 
charge and spin densities, not the orbital currents. The other is the use of an 
approximate SDFT functional (the LSDA) in these calculations. 

In practice, the resulting orbital magnetic moments strongly underestimate 
experimental values for transition metals [34], and make wrong predictions for, 
e.g., the volume collapse in lanthanides [4] and the band gap of transition- metal 
oxides [5,35]. This behaviour was interpreted as a consequence of the inapplica- 
bility of Hund's second rule to the electron gas, on which the LSDA is based. To 
account for Hund's second rule, phenomenological orbital angular momentum 
dependent corrections to the LSDA were proposed by drawing on analogies with 
angular momentum dependent terms in the multiplet splitting of atoms [3,4]. 
Significant improvement of magnetic moments and related quantities is obtained 
from this so-called orbital-polarization approach [3-5,10,11,23,34,36-39]. 

In Sec.m we explain that this empirical success can be understood if the OP 
term is not considered a better SDFT functional than the LSDA, but rather an 
approximate CDFT functional. As such, it produces orbitals that can reproduce 
charge, spin and current densities, and thus also orbital magnetic moments. 

The original proposal for systems with open shells of 4/ electrons was to add 
the OP term (today also known as Brooks term) 



to the total-energy functional of SDFT. Here is a Racah parameter defined 
in terms of Slater integrals Fk as = (5F2 + 6F4 - 91F6)/3 [40,41], and L is 
the total orbital angular momentum of the open-shell ion defined as [4,5,10] 



in terms of a sum over occupation numbers ^nimia of single-particle orbitals. 
Expressions (fT6|) and (fT7|) are motivated by noting the presence of similar 



A£;^ = — E^L^ 



(16) 



2 




(17) 
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terms in the vector model of the multiplet spHtting of atoms with open / 
shehs [4-6,40,41], but their inclusion in an SDFT calculation is entirely ad 
hoc. Nevertheless, by adopting Eq. (|16p as a starting point, simple and empir- 
ically successful correction terms to the single-particle equations of SDFT are 
obtained [3-5, 10, 11, 23, 34, 36-39]. 

The original OP idea was subsequently extended to d-electrons, for which 
Norman proposed [5] 

^^Nd' ^ -IbL^, (18) 

where B is the Racah parameter B — F2 — 5F4 [40,41]. 

Further refinements are based on a more complete treatment of atomic mul- 
tiplet splitting [6,7] and crystal-field effects [35]. For d-electrons, Norman [6] 
proposed the expression 

AE^"^ ^ 2.25nd(5 - nd)B - 1.5L{L + 1)B, (19) 

where is the occupation number of the 3d orbitals, and for /-electrons, 

AE^f = {-2nj + Unf)E^ - L{L + 1)E^, (20) 

where n/ is the occupation number of the 4/ orbitals. Shick e Gubanov [7] 
refined Norman's argument for /-electrons, and obtained 

Ai^^'=--^[L(L+l)-6.g(G2)]^;^ (21) 

where 3(62) = [6nf{7-nf)+Aln'j{7-nff-2nj{7-nff]/{2'^3^5). This expres- 
sion by construction recovers the correct splitting of the highest-spin multiplet 
of the isolated ion [6, 7]. 

Although expressions (fTOl) , (PO)) , and (PT|) give a better description of atomic 
multiplets than Eqs. (fTB|) and (HH]), applications to solids containing transition- 
metal and rare-earth atoms [4,5,10,11,23,34,36,37,39] mostly employed the 
simpler expressions (I16|) and psp . which already lead to substantial improve- 
ment of orbital magnetic moments and related quantities, as compared to un- 
corrected LSDA. In Sec. 21 we therefore mainly focus on the simple expressions 
dm) and CH]). 

We note in passing that apparently no OP term has been proposed for p 
electrons, although the physical motivation and the mathematical argument 
leading to the OP terms for d and / electrons would apply to them, too. From 
an analysis of the multiplet splitting of the p configurations [41], we find that 
the p-electron OP term should have the form 

AEP = ^ [np(6 - Up) - L{L + 1) - AS{S + 1)] F2, (22) 

which is obtained by subtracting from the energy of the terms arising from the 
p-shell their spherical average. This expression for p electrons is valid for all S, 
in contrast to expressions for d [6] and ([^0)1 and (pij) for / [6,7] electrons. 
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which hold only for the highest spin multiplet, where S — rid/ 2 and S = nf/2, 
respectively. In analogy to Eqs. (fT6|) and (fTS]) . a simpler version, expected to 
simulate the effect of (P^ . is 

AEP' = -\f2L^- (23) 

Finally, we remark that d and /-electron OP terms of similar form to ^9]) , 
(|20|) and (j2ip have been obtained also within relativistic DFT [38] . In that work, 
as well as in Refs. [34,39], it was pointed out that there should be a connection 
between the various phenomenological orbital-polarization corrections and the 
first-principles formalism of (relativistic or nonrelativistic) CDFT. In the next 
section we demonstrate this connection. 



4 Connection between the phenomenological and 
the first-principles approach 

In this section we establish a connection between the OP terms and CDFT. This 
connection is obtained in two different ways, once by analysing single-particle 
Hamiltonians (Sec. 14. and once by considering total energies (Sec. 14. 2| . The 
method of analysis is different in both cases, but the final results are the same. 

4.1 Single-particle energies 

In the absence of external magnetic fields, the Kohn-Sham single-particle Hamil- 
tonians of CDFT and DFT take the form 

= -^y' - ^ [A.c(r) • V + V • A.e(r)] + ^;(r) + vnir) + v^v) (24) 

and 

H — + v{v) + VH{r)+vU^), (25) 
respectively. There difference is thus given by the operator 

AH = H'-H= [A,,(r) • V + V • A,,(r)] + A«,e(r), (26) 

Zmc 

where Avj-d^) — v^^r) — Vxd^)- A perturbative treatment of the first term 
of this difference was suggested in Ref. [42] and applied mainly to p electrons 
in [15, 16]. In the context of magnetism, however, our main interest is in d and 
/ electrons, for which the difference between both Hamiltonians may be too 
large to justify a low-order perturbation treatment. Moreover, the functional 
used in [15, 16], as well as in many other applications of CDFT [21-23], is 
based on the electron liquid, which at stronger magnetic fields displays quantum 
oscillations [30,31], which are incorrect for finite atomic and molecular systems. 
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Independently of the size of the various terms, the identity V- [n(r)Aa;c(r)] = 
always holds [1,2], and can be used to write Ai? in spherical polar coordinates 
(r, 9, if) as 



AH = At;^,(r) 

2mc [ \ dr r 89 rsm9 dip J n{v) 



(27) 



This is the complete set of (nonrelativistic) corrections to DFT predicted by 
CDFT in the absence of external fields. 

We can now compare this to the various ad hoc corrections listed in Sec. 13.21 
According to Janak's theorem [43], the DFT Kohn-Shani eigenvalue is obtained 
from the total energy by differentiating with respect to the occupation number, 

dE 

f-nlmi<y — • (28) 

If this relation, which also holds in SDFT and CDFT, is applied to the SDFT 
functional augmented by the orbital polarization term (fTB)) . and the Racah 
parameter is kept fixed during the differentiation, it follows that the ad hoc OP 
correction to the total energy corresponds to a correction to the single-particle 
energies of the form 

Ae^ = -E^Lmu (29) 
which can in turn can be interpreted as a result of the operator 

AH^ = -E^Li,, (30) 

added to the single-particle equation of SDFT. Similarly, 

AH^'^' = -BLi,, (31) 

AHP' = -F2Lh, (32) 

^jjNd ^ 2.25 (5 - 2nd) B -'iB (^+^ 4, (33) 

AH^f = {-Anf + 14) E^ - 2E^ [L + ^ ) 4, (34) 



and 



2) " 22335 

(7 - n/) (6 - S2n]) + (7 - n/)^ (82?i/ -I- Qn)) ~ 6n} (7 - n/)^ - 6n/ E^. (35) 

We note that if the Racah term is not treated as a number, but as a matrix 
element involving single-particle radial orbitals [40,41], additional terms appear 
on the right-hand side of these equations. (We return to this issue at the end 
of Sec. 14.21 from a slightly different point of view.) 
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Since Iz — ~id/d(p, the single-particle operators resulting from the OP 
terms, Eqs. ((SO]), and can be cast in the form of Eq. ^7^, with 



A°^(r) = -— i?;rL sin (36) 

along with 

«jf^(r)=«.,(r), (37) 

and A^f'''{r) = A^f'^{r) = Vn(r)-A£f (r) = 0. In these equations ip is the unit 
vector in the azimuthal direction, and Ri = E^,B,F2, for /, d, and p electrons, 
respectively. 

This way of writing the OP correction shows that it enters the KS Hamil- 
tonian in the same way as the xc vector potential of CDFT. To complete the 
proof that it really is such a potential, we need to write it as a functional of the 
current density, A^/'[jp(r)], which is easily accomplished by means of Eqs. (|14p 
and p7)) . 

We note that v^''^^{r) — v^d^) implies Av^^{r) ~ 0, which may be inter- 
preted as a weak coupling approximation, in which the self-consistent effect of 
currents on the electric potential is considered small. The neglect of the radial 
and polar components of Aj;c(r) is correct for systems with spherical symmetry, 
where Vn(r) ■ A^c also vanishes. The orbital-polarization terms can thus be 
identified as weak-coupling, spherically symmetric approximations to the self- 
consistent xc vector potential of CDFT. 

This conclusion was anticipated in Refs. [34,39], where it was, however, ar- 
gued that the physics described by the Brooks term was distinct from that of 
CDFT. Our present point of view is different: CDFT is a general framework 
for describing all nonrelativistic effects of orbital magnetism, and to the ex- 
tent that the orbital polarization terms are justifiable, they must be specific 
approximations to the general framework of CDFT. 

The other orbital-polarization terms can also be identified as approximations 
to A^c- Eqs. (iMl) and ^ lead to 



and 

for d electrons, and 
and 



Af/(r) = -3^5 (^L + r sin ^ (38) 
v','f''{r) = vU-r) -f 2.25 (5 - 2nd) B (39) 



A^/(r) = -2^i?^ (L + \]rsine^ (40) 



vlff{v) = vU^) + {-Anj + 14) (41) 

for / electrons. Interestingly, the scalar xc potentials are also modified by 
these proposals. However, the modification is a constant shift, independent of 
position, whereas Ai;a;c(r) = v'^^r) — U2;c(r) should, generically, depend on r. 



11 



Finally, the proposal of Shick e Giibanov implies 

Aff (r) = -3^£;3 rsine ^ (42) 

and 

v^'^^^ir) = v^c{r) + — [21 + 2003n/ - 1890n^ + 670n} - 105nj + 6n}] . (43) 

All these expressions for Axd'r) and v^^(r) can be implemented self-consistently, 
by recalculating the Racah parameters and the occupation numbers at every 
iteration, although such a self-consistent recalculation is not consistent with 
treating these parameters as numbers when deducing the eigenvalue correction 
from the total-energy correction. 



4.2 Total energies 

The connection between the CDFT xc potentials and the orbital polarization 
corrections, summarized in Eqs. ([36]) and (|37p was obtained by comparing oper- 
ators in the single-particle equations of CDFT and of the SDFT -1- OP approach. 
The operator in the SDFT + OP equations was obtained by means of Janak's 
theorem, from partial differentiation of the total energy with respect to occu- 
pation numbers. 

Alternatively, one can follow the established DFT way to obtain potentials 
directly from total energies, by variational differentiation with respect to densi- 
ties, according to [1,2] 



and 



(r) 



6n{r) 



(44) 



(45) 



Application of these definitions to the DFT total energy expression, augmented 
by the OP terms HH), ([TH]) or (USD, leads to 



r = 



SExcjn] _ 1 5 
6n{r) 2 Sn{r) 

Vxc(r) - -- 



2 5n{Y) 



(46) 



and 



.OP 



(r) 



SExcjn] _ 1 
'^jp(i") 



2<5jp(r) 



1 



-If \2 

2^ <5j,(r)- 



(47) 
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These equations can be simplified by using Eq. ([T^ to calculate 5 {L z) / 5ip{Y) . 
Since within CDFT n(r) and jp(r) are independent variables, we arrive at 



and 



mc 
he 



Ri rsm9 {L,) 'p-7r{L,Y^ 



SRi 



2e 



(48) 



(49) 



If we again treat the Racah parameters as numbers, and not as self-consistent 
functionals of the densities, we finally obtain 



<f^(r)=«.,(r) 



and 



.OP 



(r) 



he 



Ri r sin 9 {Lz)ip. 



(50) 



(51) 



These are the same relations obtained in the preceding section from Janak's 
theorem and partial differentiation with respect to occupation numbers. Just as 
there, the Racah parameters must be treated as fixed numbers in order to arrive 
at the standard form of the OP term. We note that if the Racah coefficients 
are differentiated consistently also in the derivation via Janak's theorem, one 
obtains, by application of the chain rule in the form 



6Ri dn{r) , SRi ajp(r) 



5n{r) dj Sjp{r) dj 



(52) 



the same extra terms appearing on the right-hand side of Eqs. and 

To the best of our knowledge, consequences of these additional term have 

never been systematically explored. 



5 Conclusions 

The /-electron orbital polarization term proposed by Brooks at al. [3,4], its gen- 
eralization to d electrons [5], and their refinements and generalizations [6,7,35], 
have been identified as specific approximations to Aa;c(r) and v^^{r) of CDFT. 
This identification provides formal justification for their introduction into the 
Kohn-Sham Hamiltonian of (S)DFT: by using an orbital polarization correc- 
tion one actually simulates certain aspects of CDFT within the computational 
framework of (S)DFT. 

The existence of such a connection implies that the OP terms should not be 
interpreted just as corrections to the LSDA of SDFT, although this may have 
been their original motivation, but also (and perhaps principally!) as a step 
from SDFT to CDFT. Hence, from the present point of view, the OP terms 
address both issues facing the calculation of orbital moments in SDFT that 
were mentioned at the beginning of Sec. 13.21 correction of electron-gas-based 
xc functionals for Hund's second rule, and substitution of the SDFT orbitals 
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by CDFT orbitals that reproduce the orbital currents. We stress that this also 
holds for the p-electron OP terms (|22p and ((23)) . proposed in this work. 

Explicit expressions for A^c and v^^{r) that can be used in CDFT are scarce, 
and mostly based on many-body calculations of the energy of the uniform elec- 
tron gas in certain ranges of external magnetic fields [2,29-33]. Explicit, albeit 
approximate, expressions for A^c and v^^{r) that by construction go beyond 
local-density approximations may be useful in applications to finite systems, 
and constitute starting points for the development of more refined functionals. 

As an example, by combining ([5T|l with we obtain the explicit current- 
density functional 



The empirical success of the OP concept in solid-state physics suggests that 
this may be a useful and reliable expression for A^c of molecular and solid sys- 
tems containing open-shell atoms. As this expression is nonperturbative, and 
not based on the electron liquid, it is expected to be more appropriate for sys- 
tems containing open-shell atoms than the linear-response LDA [2,15,16,21-23]. 
This expression also does not suffer from quantum oscillations that stem from 
the electron liquid but are incorrect for the atomic and molecular systems 
[30,31]. It may therefore constitute a useful starting point also for applica- 
tions of CDFT to molecular systems. 
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